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Abstract. We construct a large family of exact solutions to the hyperbolic 
system of 3 equations of ideal granular hydrodynamics in several dimensions for 
arbitrary adiabatic index 7. In dependence of initial conditions these solutions 
can keep smoothness for all times or develop singularity. In particular, in the 
2D case the singularity can be formed either in a point or along a line. For 
7 = — 1 the problem is reduced to the system of two equations, related to a 
special case of the Chaplygin gas. In the ID case this system can be written 
in the Riemann invariant and can be treated in a standard way. The solution 
to the Riemann problem in this case demonstrate an unusual and complicated 
behavior. 



1. Introduction 

The motion of the dilute gas where the characteristic hydrodynamic length scale 
of the flow is sufficiently large and the viscous and heat conduction terms can be 
neglected is governed by the systems of equations of ideal granular hydrodynamics 

This system is given in R x R™, n > 1, and has the following form: 

(1) d t p + div x (pv) = 0, 

(2) d t (pv) + Div x (pv ® v) = -V x p, 

(3) d t T + {v, V X T) + (7 - l)Tdiv x v = -ApT 3 / 2 , 

where p is the gas density, v = (vi, ...,v n ) is the velocity, T is the temperature, 
p = RpT is the pressure (the constant R is a adiabatic invariant, for the sake 
of simplicity we set R = 1), and 7 is the adiabatic index, A = const > 0. We 
denote Div x and div^ the divergence of tensor and vector with respect to the space 
variables. The only difference between equations ([TjH© and the standard ideal 
gas dynamic equations (where the elastic colliding of particles is supposed) is the 
presence of the inelastic energy loss term — ApT 3 ^ 2 in ([3]). 

The granular gases are now popular subject of experimental, numerical and 
theoretical investigation (e.g. [3], [2], [E] and references therein). The Navier- 
Stokes granular hydrodynamics is the natural language for a theoretical description 
of granular macroscopic flows. A characteristic feature of time-dependent solutions 
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of the continuum equations is a formation of finite-time singularities: the density 
blowup signals the formation of close-packed clusters. 

System (HJ) — (I3J) can he written in a hyperbolic symmetric form in variables 
p, v, K = pp^ 1 and the Cauchy problem 

(p,u,Tp 1_7 )| t=0 = (p , vo,7opJ~ 7 ) 

is locally solvable in the class of smooth functions. 

System (JTJ) — (J3J) has no constant solution except the trivial one (p = 0). Another 
trivial solution is v = p = T = 0, p(t,x) = po(x). At the same time there exists a 
solution 

(4) p,v,p = const, T-T(t) = (^+T(0)- 1 /2)-2 j 

where T(0) is the initial value of temperature (the Haff's law). This solution is 
called the homogeneous cooling state. 

Here we are going to construct new exact solutions to the ideal granular hy- 
drodynamics with a concentration property and to compare them with the known 
family of solution obtained earlier in [5] . 

2. Family of exact solutions in ID [5] 

The authors employ Lagrangian coordinates and derive a broad family of ex- 
act non-stationary non-self-similar solutions. These solutions exhibit a singularity, 
where the density blowups in a finite time when starting from smooth initial con- 
ditions. Moreover, the velocity gradient also blowups while the velocity itself and 
develop a cusp discontinuity (rather then a shock) at the point of singularity. 

System (Q]) - (J3J) in the Lagrangian coordinates takes the form 

d fl\ _ dv dv_^_dp_ dp = dv 3/2 1/2 
dt \p J dm ' dt dm ' dt dm 
and under certain assumptions can be reduced to 

d 2 p 2 A 



dm*--** " = yV2- 

X 

Here m(x, t) = J p(£, t) is the Lagrangian mass coordinate. The solution is 
o 

the following: 

oiim ) 

(5) p — 2A cos(pm), A — const, p(m,t) = 1 — . 

(1 — \it^J Ap(m, 0) cos pm) 2 

The rate of concentration at the maximum point of density as t — > is p(0,t) ~ 
const(£* — i)~ 2 , the behavior of solution at different moments of time and formation 
of the singularity is presented at Fig.l. 

3. Solutions with a constraint 

Let us introduce a new dependent variable z(t,x) as follows: z = p — 4>(t)T~i , 
<j){t) is an arbitrary differentiable function. Thus, in the variables z, T, v the system 
(III)-© takes the form 

(6) d t z + div x (zv) - | 4>{t) z + ( 7 + 1) 4>{t) T-i div x v + (20' (t) + A0 2 (t)) T~5 = 0, 
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FIGURE 1 . Finite mass solution 8 



(7) 



d t v + (v,V x )v 



1 



-V x {zT + 4>{t)T*), 



(8) 



d t T + (v, W X T) + (7 - 1) T div x v = -Az T 3/2 - A <j>{t) T. 



We consider a particular class of solutions characterized by property z — 0. There 
are two possibilities: 

• 7 = — 1 and (f>(t) is a solution to ODE that can be immediately solved: 



(9) 



(10) 



v(t,x) — a(t)x + (3(t) and 4>{t) is a solution to ODE 



The first possibility is the case of Chaplygin-like [3] gas, where the state equation 
is chosen as 



(11) 



V = Po - P 



, po = const > 0. 



The system (fT1)-(f5)). (fiTj) with the constraint z = can be reduced to a couple of 
equations 



(12) 



dtv + (v, V x ) v = T* V X (T2), 



(13) 



8 t T^ + (v, V X T^) - T^div x v = ~ (j>{t) T5, 



where <p(t) is given by ©■ 

In the ID case this system as any system of two equations can be written in the 
Riemann invariants, this allows to apply the technique usual for gas dynamics, we 
will do this in Sec l3.2l 

In the second case, for an arbitrary 7 the equation ((6]) can be satisfied only 
for v(t, x) — a(t)x + f3(t), where <p(t) solves (fT0|) . This case will be considered in 
SecO 
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3.1. Solutions with uniform deformation, arbitrary 7. It is known that for 
usual gas dynamics equations the solutions with linear profile of velocity v{t,x) — 
a(t)x + Pit), where a(t) is a matrix n x n and j3{t) is an n - vector, a; is a radius- 
vector of point, constitute a very important class of solutions |10j . For the system of 
granular hydrodynamics these solutions give a possibility to construct a singularity 
arising from initial data. 

First of all from (fT2| ([T3"l) we get that in this case T has to solve the system 

(14) (d t a(t) + a 2 (t))x + (d t /3(t) + a(t)(3(t)) = ~ V X T, 

(15) d t T + ((ax + P),V x T) + { 7 -l)Ttra(t)=-A<f>(t)T, 

and the structure of the field of velocity requires a special structure of the field of 
temperature, namely, 

(16) T(t, x) = x T A{t) x + (B(t),x) + C(t). 

Thus, we get a system of 3 " + 2 5 "+ 4 nonlinear differential equations for components 
of the square matrix a(t), the square symmetric matrix A(t), vectors /3(t) and Bit), 
the scalar functions C(t) and 4>(t), namely 

(17) a'(t)+a 2 (t)+A(t) = 0, f3' (t) + 2a(t)f3(t) + ^B{t) = 0, 

(18) A'(t) + 2A(t)a(t) + (7 - l)tr a{t)A(t) + A<f>(t)A(t) = 0, 

(19) B'{t) + 2A{t)P{t) + B(t)a(t) + (7 - l)tr a(t)B(t) + A(f>(t)B(t) = 0, 

(20) C'(t) + (B(t), p{t)) + (7 - l)tr a(t)C(t) + A0(t)C(t) = 0, 

and (jlOl) . This system can be explicitly (in the simplest cases) or numerically 
integrated, one can study its qualitative behavior. The component of density can 
be found as 

(21) P[t ' X) = ix T Ait)x + iBit),x)+Cit)) 1 /^ 

pit,x) ~ const * \x — xo| _1 in the point Xq of the singularity formation. Therefore 
the singularity is integrable for n > 1. Nevertheless, the total mass is infinite for 
this solution, since J pdx diverges as \x\ — > 00. 

Let us consider the simplest non-rotational case: A(t) = a(t)I, a(t) = a>i(t)I 
Bit) — 0, Pit) = 0, where I is the unit matrix. The system above comes to 4 
equations: 

n A 

(22) <t,>(t) + -i 7 +l)4>it) ai it)--<p 2 it), ai(t) + a?(t)+o(t)=0, 

(23) a 1 it) + ((2 + n(7 - l))o x (t) + A0(t))o(t) = 0, 

(24) C"(t) + (n(7 - l)ai(t) + A^(t))C(t) = 0. 

We are going to find asymptotics of the solution at the point t — t* > of the 
singularity appearance. 

Systems (|2"2")) - (|2"41 is a polynomial system 

(25) x = f(x), f:R fc ^R fc , 
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and we can study the occurrence of blow-up analyzing the solutions locally around 
their movable singularities using a set of methods based on the construction of local 
series. Following [5J, [BJ, we build local series (<J/-series) of the form: 

x = $(\,s,t) = Xt s (1 + /i(r,lnr)), 

where r = t» — t and h(r; lnr) is a power series in its argument which vanishes as 
t — > 0. The notation Xt s refers to the vector whose i-th component is XiT Si . In 
order to obtain the leading behavior Xt 3 of the solution around i* we look for all 
negatively quasihomogeneous truncations f of the vector field f = f + f such that the 
dominant behavior x = Ar s , A G C fe is an exact solution of the truncated system 
x = f(x) and 

f (Ar 5 ) ~ AV+ 1 " 1 , seQ k , §i>0, 
as t — > 0. Each truncation defines a dominant balance (A, s) and every balance 
corresponds to the first term Xt s in an expansion around movable singularities. 
For such an expansion to describe a general solution, the Vf-series must contain 
k — 1 arbitrary constants in addition to the arbitrary parameter t* . The position in 
the power series where these arbitrary constants appear is given by the resonances. 
They are given by the eigenvalues of the matrix R: 

R = -.Df(A) - diag(s), 

where Di(X) is the Jacobian matrix evaluated on A. The resonances are labeled 
7"j, i = 1, k with ri = — 1. Each balance defines a new set of resonances. 

Theorem 3.1. (]6j, T.3.8) Consider a real analytic system h25\) and assume that 
it has a balance (A, s) such that rj > for all j > k — m + 1, 1 < m < k and 
X £ 1". Then there exists a m- dimensional manifold S™ C M. k of initial conditions 
leading to a finite time blow-up, for all Xq € S™ , i.e. there exists € M+ for which 
|x(£;xo)| — > oo as t — > 4*. 

Corollary 3.1. If n > 2 anc! 7 > — 1 + — , £/iert £/iere exists an open set SI 0/ 
initial data x = (4>(0), Qfi(0), a(0), C(0)) smc/i i/ia< i/ie components <fi(t) and cti(t) 
of solution to the system ^28 \) - (E$\) blow up within a finite time for all Xq € 17. 

Proof. To find main terms of asymptotic at the point of singularity we consider a 
negatively quasihomogeneous truncation of the system (|2^t - (P^|) . namely 

(26) ^(t) = -^(7+l)^(*)oi(t)-^(t), 

(27) a[(t) = -af(t), 

(28) a'(i) = -(2 + n(7 - l))ai(t)a(t) - A<f>(t)a(i), 

(29) C"(t) = -n(7 - l)ai(i)C(t) - A<p(t)C(t). 
The solution to the above system is the following: 

(30) <j>(t)=- nin+ ^~ 2 (t*-t)-\ a 1 (t) = -(U-t)-\ 

A(t) = A (t* -t) 2(n - 2) , G(i) =C (t*-i) 2( ™- 1) , A , C = const., 
s = rfza.g(-l, -l,2(n - 2), 2(n - 1)), A = (-1, + 1 ~ 2) A 0) C ). 
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The resonances, computed for this balance are (71(7+1)— 2, —1, 0, 0). Theorem l3.ll 
result that there exists a manyfold Sq £ R" such that for c/>(0), ai(0) the respective 
solution to (j2"2")l blows up and has (|5T?|) as a main term of asymptotics. Other 
components of solution to (f22j) - ([24]) can be found from linear with respect to a(t) 
and C{t) equations (|28|) and (|29l) for any initial data (it makes sense to consider 
o(0) > 0,(7(0) > 0). □ 

Remark 3.1. The rate of growth of the maximum of the density as t — > t» is 
p(t,0) ~ constat-**) - ". 

Remark 3.2. If i* < 0, then an analogous consideration shows that there exists an 
open set of initial data such that the solution to system (|22 p ~ (|24p remains bounded 
for all t > 0. 

Theorem 3.2. It n = 1, then for any t* > i/iere exists a family of exact solutions 
to system I22fy - \21$ depending on parameters (ao, Co), blowing up as t — » t». 77ws 

family is physically reasonable for «o € (— 1, ; ^j), Co > 0,7 > 1. for these 

solutions the the maximum of density has the asymptotics p(t,0) ~ const(i — £*) Q ° 
as t — > t* . 

Proof. It can be readily checked that the balance s = (— 1, — 1, — 2, — 2(ao + 1)), 
A = (— (2+(7+l)ao)/A, a , — ao(ao+l), Co) gives an exact solution. The restriction 
on the parameters ao and Co follows from the positivity of the expression under 
the square root in (|2ip and the positivity of </>(0). □ 

Remark 3.3. The maximum of density as t — > t* grows slower than for the solution 
obtained in [5]. 

Fig. 2a presents the results of numerical computations in 2D based on system 
([17)1 The initial density has the form ([2"T]l. Fig. 2b shows the density near the 

blow-up time for an — ai2 < 0, a\i = C121 = 0, with a concentration in a point. 
Fig. 2c shows the density near the blow-up time for an < 0, «22 = c*i2 = 0^21 = 
with a concentration along a line. The computations demonstrate a complicated 
behavior of solution. In particular, a vorticity can prevent the singularity formation. 




Figure 2. a,b,c 
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3.2. Chaplygin gas, n = 1. The model of gas dynamics with the pressure given 
by (fill) i s known as the Chaplygin gas. The Chaplygin gas is now considered as 
a possible model for dark matter-energy (T). This system can also be seen as the 
one-dimensional version of the Born-Infeld system, a non linear modification of the 
Maxwell equations, designed by Born and Infeld in 1934 to solve the electrostatic 
divergence generated by point particles in classical electrodynamics. The Chaplygin 
system is known to be hyperbolic, linearly degenerate, weakly stable [12] . Recently 
this system attracted a lot of attention, e.g [2J, [TTj . 
The system p2[) . (IT31 can be reduced to 

(31) d t p + div x (pv) =0, d t (pv)+mv x (pv®v- =0, 

P 

recall that p = 4>{t)T~^. It is similar to the Chaplygin gas system, the only 
difference is in the known multiplier (j)(t), for the Chaplygin gas <j) — const. For 
n = 1 (1311) can be written in the Riemann invariants as 



(32) d tS + r d x s = ^±(r- S ), d t r + sd x r = ^-{s - r), 

where s — v~T^ , r = v + T^ , <j>(t) is given by ©■ This system is linear degenerate, 
therefore provided the solution is bounded there is no possibility for the gradient 
catastrophe. 

Theorem 3.3. The solution to the Riemann problem for \31}) ( \3%fy ) in ID with 
data 



(33) (v,T) 



(vl,T l ), x<0, 
(vr,T r ), x>0, 



in the case 

(34) v L > v R 
contains a S-singularity in the component of density. If 

(35) v L >v R + Tl+Tl, 

then the S-singularity formation begins from the initial moment of time. 

Proof. Since the system is linear degenerate, the jumps are contact discontinuities 
and move along characteristics. The solution is based on the cooling state ([4]). If 

(36) v L <v R + Tl+Tl, 
then the solution is 

(v L ,T L (t)), x<x-(t), 

(37) (v,T) = { (v M (t),T M (t)), x-(t) <x<x+{t), 

(v R ,T L (t)), x>x+(t), 



with c = r^O), T L {t) = jh^, T R (t) 



c'T R 



(ft+c) 3 ' (#*+c) 2 ' 



2cT 2 A 2cT 2 A 

x-(t) - v L t - —± ln(-t + 1) x+ {t) = v R t + ln(— t + 1), 

M VL + VR + cjTl-Tl^t + c)- 1 

VM(t) - , 
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1 M\ l ) — 2 

If vl < vr, then X-(t) < x+(t) for all t > and the solution to the Riemann 
problem is given by (|37p . If Vl < vr, then there exists a moment f* > such that 
ir_(t*) = x+(t*). Moreover, in the moment t** > the component Tm vanishes 
and i** < t*. Thus, we have to construct a new solution starting from t**. 

We are going to introduce a (5-singularity in the density concentrated on the jump 
analogously to [2] . To find a <5-type singularity solution we have to use the system 
in its conservative form pip . Let us denote x*(t) the position of the singularity 
and look for a solution in the form: 

(38) p(t, x) = p- + [p]H(x - x*(t)) + 6(t)6(x - x*(t)), 

(39) p(t,x)v(t,x) = pv- + [pv]H(x- x*(t)) + ip(t)6(x - x*(t)), 

(40) p{t, x)v 2 (t, x) = pv 2 _ - +[pv 2 ]H(x - x*(t)) + V(t)6(x - x*(t)), 

(41) T(t>x)=T- + [T]H(x-x*(t)), r^p- 1 , 

[/] = /_|_ — /_ , /_)_ and /_ are the limits of an arbitrary function from the right and 
from the left side of x*(£), respectively, H is the Heaviside function. From (j3"Tj) we 
get 

(42) fcW -H!£M 

[p\ 



(43) e(t) = M[pv] 2 -[p][pv 2 ])t\ + 



fi = t — to, to and xo are the moment and the coordinate of the 5-singularity 
formation. Expanding the expression under the square root at t = to we find the 
necessary condition for the beginning of the concentration processes: 

(44) H 2 > [T 1//2 (t )] 2 - 

If the initial data are such that the inequality opposite to (|36[) holds, i.e. vl > 
vr + Tl /2 (0) + T^ /2 (0), then v L > v R and 



(v L - v R ) 2 > (T^ /2 (0) + < 2 (0)) 2 > (T^ /2 (0) - T^ /2 (0)) 2 



such that the condition (|44j) is satisfied. Since in this case the solution consisting 
of two contact discontinuities is impossible, the only reasonable solution is given by 

dssh dm, dm. 

If ul > ur, we can define again the solution of form (I3"51) - (|4"Tj) . (|4"2")l . In principle, 
initially the condition (|44[) may fail and this solution can exist beginning from some 
t, < t < nevertheless, this moment t always exists. Further, at least for t>t* 
this solution this solution is stable. Let us extend this solution back < t < i*. 
This means the assumption that the segment [x_(f),x+(t)] shrinks into the point 
x*(t) at the moment t**. It can be shown that the velocity of the singular front 
xjt) -> ! H Ea as i -> cx). □ 

Remark 3.4. The question on uniqueness of the solution after the "shrinking" is 
open. 
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Remark 3.5. In [9] for any spatial dimensions a simple family of solutions to the 
system ([I]) - ((3|) having a singularity in the density whereas other components are 
continuous is constructed. Moreover, a family of self-similar solutions in ID was 
found. 
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